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■ Abstract 



A generalization of the Gram-Schmidt procedure is achieved by providing equations 
for updating and downdating obhque projectors. The work is motivated by the problem 
of adaptive signal representation outside the orthogonal basis setting. The proposed 
I techniques are shown to be relevant to the problem of discriminating signals produced by 

' different phenomena when the order of the signal model needs to be adjusted. 



■ 1 Introduction 

> . 



An operator (9 is a projector if it is idempotent, i.e., if it is endowed with property O = O. 



O ' The projection is along (or parallel to) its null space M{0) and onto its range TZ{0). This 
5 ■ entails that Ov = v ioi v e 7^(C») and Ov = Q for v G U{0). If the subspaces 7^(C) and M{0) 
\Q . are orthogonal the operator is called an orthogonal projector, which is the case if and only if 
^ i C is self-adjoint. Otherwise it is called oblique projector. 

i-g \ Oblique projectors, though introduced early [1,2], have received less attention than orthog- 

■ onal projectors. Nevertheless, quite recently there has been a renewed interest in relation to 
their properties and applications [3-12]. In particular, oblique projectors have been shown to 
be of significant relevance to signal processing techniques [13-16]. The present effort is very 
^ ! much motivated by problems arising in the area of signal representation outside the traditional 
^ \ orthogonal basis setting [17-20]. In such a context a signal /, represented mathematically as 
an element of a vector space, is approximated as a linear expansion of the form 

k 
i=l 

The vectors Vi in ([1]) are sometimes sequentially fed or chosen according to some optimality 
criterion. In such situations one needs to be in a position to effectively adapt the coefficients 
of the linear superposition so as to account for the possibility of changes in the model. This 
may entail a) increasing the order k of the model by incorporating new terms in the expansion 
b)reducing the order by eliminating some terms in the expansion c)replacing some of the vectors 
in ([1]) by different ones. 

Assuming that the signal space is an inner product space, for given in ([1]) to be the 
best approximation of a signal / in a minimum distance sense, the coefficients in ([1]) should 
be calculated in such a way that /'^ is the orthogonal projection of / onto span{fj}^^^. This 
is the main reason for the popularity of orthogonal projectors in the context of approximation 



1 



techniques. Nevertheless, suppose that the observed signal is produced by the interference of 
two phenomena so that the model ([1]) becomes 

k n 

'^CiVi + ^^diWi. (2) 

1=1 1=1 

If one were interested in discriminating the phenomena by splitting the signal, the component 
in span{fj}^^^ could be obtained by an oblique projection operation mapping the other com- 
ponent to zero. There is a broad range of applications in which this procedure happens to 
be of assistance [13]. Thus, we felt motivated to find recursive equations for adapting oblique 
projectors. Some of the equations to be proposed here are inspired by our previous work on re- 
cursive biorthogonalization for orthogonal projectors representation [21,22]. We have recently 
been made aware that such a work is closely related to earlier one on recursive generalized 
inverses [23-26]. 

In spite of the fact that for most numerical implementations a projector is represented by 
a matrix, we prefer to think of projectors as operators acting by performing inner products. 
An important reason for this choice is the following: The equations can thereby be applied in 
general inner product spaces and comprise two very important cases in particular. Namely, the 
Euclidean inner product space, where a projector is indeed a matrix, and the space of functions 
of finite 2-norm. We like to see the proposed recursive equations as generalized Gram-Schmidt 
like procedures for generating sequences in inner product spaces. Such sequences give rise to 
oblique projectors onto nested subspaces and, of course, to orthogonal projectors as special 
case. 

The paper is organized as follows: Section [2] introduces the notation along with a discussion 
on the general construction of oblique projectors. Section |3] provides the recursive equations 
for stepwise updating/downdating of such projectors. Applications are illustrated in Section H] 
by i)recovering a simulated X-ray diffraction peak from a background and ii)filtering impulsive 
noise from the register of the motion of a system consisting of the superposition of damped 
harmonic oscillators. The conclusions are drawn in Section [51 



2 Oblique projectors 

As already mentioned we will work in a general inner product space 7i, where the square norm 
is induced by the inner product that we represent as (■,■). Given two closed subspaces, 
V ETi and W"*" G Ti, such that Ti = V + W"*" and Vfl W"*" = {0}, the oblique projector opersitoY 
onto V along W"*" will be represented as -Eyw^ ■ Then -Eyw^ satisfies: 

-Eyw-Lf = V, for any t> G V 
Eyy^^w = 0, for any w G W"*". 

In the particular situation in which W"*" happens to be the orthogonal complement of V in 7i, 
i.e. if W"*" = V^, the operator is self-adjoint and represents an orthogonal projection onto V. 
We emphasize this special case by using the particular notation i^yv^ = Pv- In the sequel the 
orthogonal projector operator onto a subspace, say the subspace X, will be indicated as Px- 

Let us assume that in general V = spanjwjj^L]^ and W = spanjuj}^^^, with W the orthogonal 
complement of W^. Denoting as Cj, z = 1, . . . , the standard orthonormal basis in C*', i.e.. 
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the inner product (ej, Sj) = 6ij with 6ij equal one if i = j and zero otherwise, we define the 
operators V and W as 



k 



^ = Zl^^^^^'')' W = J2'^i(^i^■)■ 



^=l 1=1 

Thus the corresponding adjoint operators W* and V* are 

k k 

1=1 i=l 

The operations {vi, ■) and {ui, ■) indicate that V* and W* act by performing inner products in 
Ti. The inner product is defined in such a way that for f E Ti and c a complex constant the 
mapping V*cf produces a vector in of the form V*cf = cY^\=i^i{vi, f) . The operation 
(cj, ■) indicates the inner product in C^, thereby for r G C'^ the mapping Vcr yields a vector in 

V of the form Vcr = c Yli=i '^ii^i: Note that the matrix representation of W*V has elements 
given by the inner products {ui, Vj), i, j = 1, . . . , k. The operator 

where (■)^ denotes the Moore-Penrose pseudo-inverse, is known to be the oblique projector onto 

V along W"*" [15]. The particular choice 

Ui = Vi- Py^^Vi = PwVi, i= l,...,k (3) 

produces the expression for Eyy^± used in signal processing applications [13]. Certainly, setting 
W = PyjV one has the convenient equation 

^vw^ = Viy*P^V)^V*P^ (4) 

that we adopt hereafter. 

Amongst the many properties of oblique projectors that have been studied we shall recall 
only the basic property needed for our purpose. It follows by applying Pyy on both sides of (jl]), 
i.e., 

Pw^vw^ = PwV{y*PwVYV*P^. (5) 

Since 1/*P^ = V*Pw and V*PwV = V*P^V, with the substitution A = Py^V the right hand 
side of ([5]) takes the form y4(y4*y4)M*. Such an operator is the orthogonal projector onto 1Z{A). 
Consequently, 

PwEvw^ = Av, with W = n{PwV) = span{P^ui}ti- (6) 
By denoting Ui = Yl'j=i dlj'^j with gjj the element of matrix {V*Py^;Vy , we can express 

Eyy\;± as 

k 

E^^± =^Vi{ui,-). (7) 

i=l 

Furthermore, from ([6]), ([7]), and ([3]) 

k 

Pyv = PwEvy^i- = ^ Ui{ui, ■). (8) 

i=l 



Because Pw is self-adjoint span{Mj}f^^ = span{Mj}*Lj^ = W, and vice versa. On comparing ([7]) 
and ([8]) we see that the dual vectors Ui are the same. This is of enormous assistance to derive 
the equations for adapting oblique projectors so as to account for the updating or downdating 
of the projecting subspace V. This will allow us to give the proofs of the proposed recursive 
equations either by verification or by induction. 

Remark 1. It is appropriate to stress at this point that if we chose W"*" = V"*" we would have 
Ui = Vi, i = l,...,k and consequently spanjujj^^^^ = span{wj}jL]^. Hence for such special 
situation -Eyv-L = Pw = Pv Oi'nd all the recursive equations of the subsequent sections would 
give rise to orthogonal projectors. 



3 Constructing recursive equations 

In this section we provide the equations for updating and downdating oblique projectors in 
order to account for the following situations: 

Let us consider that the oblique projector -Ey^w-i- onto the subspace Vk = span{t>j}f^]^ along 
a given subspace W"'" is known. If the subspace Vk is enlarged to V^+i by the inclusion of 
one element, i.e., V^+i = span{t>j}(^^]^, we wish to construct E'y^^-^vv^ from the availabihty of 
-^VfeVV-L- On the other hand, if the subspace Vk = span{t>j}f^^ is reduced by the elimination 
of one element, say the j-th one, we wish to construct the corresponding oblique projector 
-^Vfc\j W-i- from the knowledge of Ey^-^^±. The subspace W"'" is assumed to be fixed. Its orthogonal 
complement Wfe in TCk = Vk + W"*" changes with the index k to satisfy TCk = Wk © VV"*", where 
© denotes an orthogonal sum whilst the former is a direct sum, i.e., Vk H W"*" = {0}. 

3.1 Updating the oblique projector Ey^y^± to E'y^^^yyj. 
We assume that E^^y^i^ is known and write it in the explicit form 

k 

^v,w^ =5Z^i(S-,-). (9) 

i=l 

Our aim is to find the vector u^k^i, and to change the vectors , i = 1, . . . , A; to u^^^, i = 1, . . . , fc, 
so as to obtain 

fc+i 

E^,^,w^ = Y.v,{v^^\-). (10) 

i=l 

We will show that the duals u^^^., i = 1, . . . , + 1 can be constructed inductively from the dual 
of a single vector. 

Lemma 1. For u\ = yi^^fp' ^^^^ ~ Py\>Vi, operator Vi{u\, ■) is the oblique projector onto the 
span of the single vector Vi along W"*" . 

Proof. From the definition of Ui (C.f. eq. ([3])) it follows that the operator Vi{u\, ■) = fi(pYj]2 > ') 
maps every vector in W"*" to the zero vector. Suppose that / is in the span of vi. Then / = cvi 
for some constant c. Since {ui,ui) = {ui,vi — Py^^vi) = {ui,vi) we have 



vi{u^, of) = CVi ' = CVi = f. 
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Moreover vi{u\, vi) {u\, ■) = vi{u\,-), which concludes the proof that vi{u\,-) is the obhque 
projector onto the span of Vi along W"'". □ 

In order to inductively construct from u\ = j^^f^ the duals u'l'^'^, i = 1, . . . , k + 1 we have 
to discriminate two possibilities 

i) Vfc+i = span{fi}('+^^ = spaii{vi}i^^ = Vk, i.e., Vk+i e Vk- 

ii) Vfc+i = span{wi}f+/ D span{vi}'i^^ = Vk, i.e. Vk+i ^ Vk- 

Let us consider first the case i). Clearly if Vk+i G Vk the corresponding Uk+i = Vk+i — Pw±Vk+i 
belongs to Wk = span{Mj}*Lj^, because Vk+i = X]i=i yields Uk+i = Yli=i'^i'^i- The proposi- 
tion below prescribes how to modify the corresponding dual vectors in order to guarantee that 

Proposition 1. Let Vk+i G Vk and vectors in ^ he given. For an arbitrary vector Uk+i G 7i 
the dual vectors "U^"^^ computed as 

u^+' = u'l-{uk+,,^)yk+i (11) 

for i = 1, . . . , and £t^} = yk+i produce the identical oblique projector as the dual vectors 
'^iy i = ■ ■ ■ , k. 

Proof. We use ( TTT]) to explicitly express -E'v^^^^i 

fc+l k k 

^Vi{u-+^,-) = ^Vi{Ui,-,) -'^Vi{u-,Uk+i){yk+u-) +Vk+i{yk+i,-) 

i=l 1=1 1=1 

= ^VkW^ - EvkW^'^k+i{yk+i, ■) + Vk+i{yk+i, ■) 

= Ey^y^± -Vk+i{yk+i,-) +Vk+i{yk+i,-), (12) 

where the last equality holds because Ey^y^^Py^i- = and Ey^y^±Vk+i = Vk+i for Vk+i G Vk- 
Hence, the left hand side of (IT^ equals Ey^-^±. □ 

The next proposition considers the case ii) 

Proposition 2. Let vector Vk+i ^ Vk and vectors u'^,i = l,...,k in be given. Thus the 
dual vectors m^^"^ computed as 

u1^' = u1-ulXl{uk+i,u1), (13) 
whereu^'^l = withqk+i = Wfe+i — Av^'^fc+i; provide us with the oblique projector Ey^^^y^;± . 

Proof. In order to organize the proof let us establish the following relations: 

{qk+i,Vi) = 0, for i = l,...,k (14) 

{qk+i,Vk+i) = {uk+i,Vk+i) - {vk+i,Pn>^Vk+i) = Wqk+iW^- (15) 

The first relation follows from the definition of qk+i and the fact that Pwk'^i = Ui for i = 1, . . . , k 

{qk+l,Vi) = {Uk+l,Vi) - {PyVk'^k+l,Vi) = {Uk+i,Ui) + {Uk+i, Pyv^Vi) " (Wfc+l, Av^^i) = 0. 
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On the other hand 
Furthermore 

Ikfc+llP = {Qk+l,Uk+l) — (qk+l, PWk''^k+l) = {Qk+l,Uk+l) 

= {(lk+i,Vk+i) — {qk+1, Pw-^Vk+i) = {qk+i,Vk+i) 

= {Uk+l, Vk+l) - {Vk+l, Avfc^fe+l)- 

We are now in a position to start the proof of the proposition by induction. From Lemma [T] we 
know that •)/||'u];|p is the obhque projector onto Vi along W-*-. Assuming that -Ey^wx 

is the obhque projector onto Vk along W"*" we will prove that Ey^^^yyi^ is the oblique projector 
onto Vfc+i along W"'". For this we need to prove that the recursive equation ( |T3l) yields the 
operator Ey^^^^y^i^ satisfying: 

ii) ^Vfe+iW^w = V, for any v E Vk+i 



iii) i^Vfc+iW^""^ = 0; fo^^ aiiy ^ ^ W . 
We begin by using f|T3l) to express ^v^+iW^ 



k+l 



1=1 i=l i=l 

For all w in W"*" it holds that Ey^yv±w = and {qk+i,w) = 0. Then from (fT6|) we conclude 
that condition iii) is satisfied. Every v G Vk+i can be written as f = Yli=i '^i'^i- Thus, from 
(fT6|) and using relations (|T^ and (|T5l) 

fc 



-^Vfc+iVV^"^ = ^ + Cfc+l-Ev^yy;Xffc+i — Cfc+l-E'y^yi;±'Ufc+i + Ck+l^k+l 

i=l 
k 

= ^ CiVi + Ck+lVk+1 = V, 



i=l 



which demonstrates condition ii). Finally, since from (1161) and f[T^ it follows that E^^^_^y^±Ey^ 



-^VfcW-Lj we have 



= -^VfeW-L - -^VfcW^^fc+i( II '^^'^^112 ' ■) + ^fc+i( II '^^'^^2 ' ■) = -^Vfc+iW-L- (17) 

□ 
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Property 1. If vectors {vi}^^l are linearly independent they are also biorthogonal to the dual 
vectors arising inductively from the recursive equation ( IT3|) . 

The proof of this property is given in Appendix A. 

Remark 2. If vectors {vi}^^^ are not linearly independent the oblique projector Ey^y^± is not 
unique. Indeed, if {11^}^^-^ are dual vectors giving rise to Ey^y^± then one can construct infinitely 
many duals as: 

k 

yi=Ui +yi-Y^ yj{vj, Ui) i = 1, . . . , k, (18) 

i=i 

where yi, i = 1, . . . , k are arbitrary vectors in Ti. 
Proof. We use ( IT8|) to write 

^Vi{yi,-) = Ey^w± + ^Vi{yi,-) -^Vi^{u'^,Vj){yj,-) 

i=l i=l i=l j=l 

k k k 

= Ey^w^ + ^Vi{yi,-) - ^^Vi{u\,Vj){yj,-) 

i=l j=l 1=1 

k k 

= Ey^w± + Vi{yi, ■) - ^ Vj{yj, ■) 
i=i j=i 

= E^,w^. (19) 

□ 

It follows from Property [1] that if vectors {f are linearly independent equation f|T8l) 
yields the unique duals yi = u^, i = 1, . . . , k. 

3.2 Downdating the oblique projector £^Vfe>v^ to £'y^ yy± 

Suppose that by the elimination of the element j the subspace Vk is reduced to Vk\j = 
span{t>j}^=i . In order to give the equations for adapting the corresponding dual vectors gener- 

ating the oblique projector Ey^. .u;± we need to consider two situations: 

i) Vk\j = span{t;j}^=i = span{fi}f^i = Vk i.e., Vj e Vk\j. 

ii) Vk\j = span{t;j}li C span{vi}i^^ = Vk, i.e., vj ^ Vk\j. 
The next proposition addresses i). 

Proposition 3. Let Ey^^± be given by ([9]) and let us assume that removing vector vj from the 
spanning set of Vk leaves the identical subspace, i.e., vj G Vk\j. Hence, if the remaining dual 
vectors are modified as follows: 



1 - (Uj, S^) ' 



uT"' = + ; '\:'Jk, ^ (20) 



the corresponding oblique projector does not change, i.e. Ey^^^y^!± = -Ey^w-i 



Proof. Let us recall that Vj G Vk\j implies Uj G Vk\j- Hence {uj,u^) ^ 1, as it is seen from 

the fact that {uj,u';) = Y!l=i{uj,u'l){ui,u)) = Eti(«i'^^)(^^«i) = Eti K^i- ^J") P. which 
implies {uj, u^) = 1 if and only if for i = 1, . . . , A; it holds that {uj,u^) = 6ij. This is not true 
if Vj G Vk\j so that we can use (l20l) to express E-];^^.^± as 

k k k I ~ k \ I ~ k \ 

^^("^ '■) = 2^<^ir) + l^ l-lv^v) 
i=i i=i i=i 



± - M • , ■ + — \ Vj {Uj , Uj)- r 



= -^VfeW-L- (21) 

□ 

Finally, proposition 4 addresses ii). 

Proposition 4. Let Ey^y^± be given by and let us assume that the vector Vj to be removed 
from the spanning set ofVk is not in Vk\j- In order to produce the oblique projector Ey^^^^± the 
appropriate modification of the dual vectors can be achieved by means of the following equation 

«i =«^-^j^F^- (22) 



v^{ulu^j){u^,-) 



[22l) we write: 








k 


k 






-E 


i=l 


i=l 


1=1 




= -^VfcW-L - Vj { 





\u,\ 



We notice that l\fi^2 is the orthogonal projector onto the span of the single vector Uj and 

denote it as P^k = -j—j^. Thus the orthogonal projector onto yVk\j can be expressed as 

j 1 1 



P'Wk\j = -Pwfe — ■ Applying this operator on the right hand side of fl23l) we obtain: 

^■^ 3 3 3 

and, since P-w^E^^y^i^ = Py^^ and P^kPy^^ = P^*, using the fact that P^kE^^y^i^ = P^kPy^^E^^y^i^ 



P,-,k it follows that 



(Av^ - P^k){Ey^y^l^ - E^^y^^P^k) - P^^ - Pr^k - Pr^k + P^ 



k 

J ' "] "J "j 

Avfe - Puk = Pwu\j = PyVk\j^Vk\jW^- 
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From the last equation and (!23l) we gather that 

k 

^w,,, ^^(^?^'' ■) - ^w,,,^v,,,w^ = 0. (24) 

i = l 

For every vector / G ?i we therefore have 

k 

Pv^.,AY.<^'^' ^ ■) - ^v,,,w-)/ = Pw,,ADf = 0, (25) 

i=l 

with AZ) = ^^=1 Vi{uf^\-) — Ey yy±. This imphes that either AD is the zero operator or 

AD/ G W"*" for every / G 7i. The latter cannot be true because from the definition of AD it 
is seen that ADf G Vk\j and by hypothesis Vk\j H W"*" = {0}. Hence AD should be the zero 
operator, which leads to the conclusion that 

k 

^v.,,w^ = E^^(^'^''-)- (26) 

i = l 

□ 

Remark 3. The case of replacing a vector in Vk, say vj by Vj, is actually equivalent to aug- 
menting the subspace Vk\j to Vk\j + Vj after the vector vj was deleted. Some implementation 
issues arise, though. In order to modify the duals as prescribed in (fT3|l we need to compute a 
vector qk\j = u'j — Py^^^M'y For the sequential enlargement of the projecting subspace, discussed 

in section \3.1\ the projector Py^^ can be sequentially constructed by means of the orthonormal 
vectors qn/\\(ln\\, n = l,...,k. Nevertheless, when replacing vectors sequentially we need to 
allow for the recalculation of the corresponding projectors. One possibility that could be consid- 
ered is the recalculation of the orthogonal vectors qn [27, 28]. An alternative approach implies to 
use of the dual corresponding to the deleted vector for orthogonalization purposes. A discussion 
concerning the implementation of such a procedure is given in [29, 30]. 



4 Applications to signals discrimination 

The examples presented in this section aim at illustrating the application of our recursive 
construction of oblique projectors for signals in L2[a , 6], the space of square integrable functions 
on [a, h\. For / and g in L2[a , 6], we define the inner product, according to the previously 
adopted convention, as 

(/,^) = / f*{x)gix)dx, 

J a 

where f*{x) indicates the complex conjugate of /. In the examples bellow all the integrals are 
numerically calculated. 
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4.1 Extraction of a X-ray diffraction peak from a dispersive back- 
ground 

Here the signal is simulated by emulating a crystallographic problem. It is assumed to be the 
X-ray diffraction intensity produced by a powder sample of a clay mineral consisting of very 
flat crystals. Each such crystal is formed by the stacking of n layers producing a diffraction 
intensity as given by [31,32] 

^ , , sin^ nx , 
In{x) = — (27) 
sm X 

where the variable x, given in radians, is related to the diffraction angle 6 according to the 
equation 

x = 27r^sin^. (28) 
A 

The parameter d in (l28|l is the effective distance between two consecutive layers and character- 
izes the material. The parameter A represents the wavelength of the incident radiation. 
We denote the diffraction intensity produced by the whole sample as /i(x). Thus, 

k . 2 

n / N sm nx , , 

/i(x)=5^c„— (29) 
sm X 

n=l 

As already mentioned n indicates the possible number of layers forming a single crystal in the 
sample. The coefficients Cn account for the proportions of crystals consisting of n layers. Here, 
for simulating the signal, the coefficients were considered to be 

c„ = + o.2e-°-i("-35)^ n = 1 . . . , 60. 

The diffraction figure /i emerges from a background that is modelled as 

3 

/2(x) = 50j]je-^(^-f). (30) 
j=i 

The combined phenomenon gives rise to the signal / = /i + /2 plotted in the left graph of 
Figure 1 on the interval relevant to the diffraction model, namely [| , 

We are interested in extracting the diffraction peak fi from the background. For this we 
will construct sequentially oblique projectors onto subspaces Vk given as 

Vk = span{— ^ — , n = 1, . . . , fc}, X G [- , — ]. 

sm X II 

The final /c- value is to be adjusted. The subspace W"*" is here 

= span{e-^(^-i), j = 1, . . . , 3}, x G [| , y ]. 

Since the order k of the diffraction model is assumed unknown, it was adjusted as follows: firstly 
the order model was sequentially increased (up to /c = 200) and then sequentially downdated. 
It was observed that the recovering of the signal was not very sensitive to the model order. 
In a range from = 50 to = 200 the approximations were totally equivalent. From = 40 
to /c = 50 changes in the approximations were noticed but the approximations could still be 
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considered 'practically' equivalent. The recovered peak /i = E^^y^±f (for k = 50) is depicted 
in the right graph of Figure 1. It happens to coincide, in the scale of the figure, with the graph 
of the theoretical one (C.f. (|29|)). 

The convenience of the proposed adaptive technique in the determination of the order of 
the diffraction intensity model is clear: otherwise when changing the k value as described above 
the whole projector would have to be recalculated for each different value of k. However, the 
advantage of the proposed technique is even more significant when, as is the case in the next 
example, overestimation of the order in the signal model may result in the failure to discriminate 
the signals. 




Figure 1: Graphs on the left: Diffraction intensity vs the variable x measured in radians and 
related to the diffraction angle through (l28l) . Graph on the right: Diffraction intensity extracted 
from the background. 



4.2 Elimination of impulsive noise 

In this case the signal is considered to be the register of the motion of a system consisting of 
uncoupled damped harmonic oscillators. The n-th oscillator is characterized by a frequency of 
I Hz. The corresponding equation for its motion x„(t) as a function of time is given as 

Xn(t) = e~'' cos^Tint), n = l,...,k. (31) 

The distribution of frequencies is considered to be = (1 + 0.7(n — 75)^)^^ so that the motion 
of the system is registered by the signal 

n=l ^ ' 

This signal (shown in the right graphs of Figure 2) is corrupted by impulsive noise, which 
represents a type of electrical noise appearing in some practical situations. The possible pulses 
are taking from the set of 400 Gaussian sparks pj(t) = e~io°o°o(*~o-°o^^-^)^, j = 1, . . . , 400. Hence 
the corresponding subspace is 

= span{e-~(*-°-°°25i)^, ^- = 1, . . . ,400}, t G [0, 1]. 
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First a random superposition of 50 pulses is added to the signal /i to simulate the noisy one 
plotted in the top left graph of Figure 2. The result after the oblique projections along the 
subspace given above is shown in the top right graph (it coincides with the theoretical 
signal /i given in fl32l) ). The left graph at the bottom of Figure 2 corresponds to a different 
realization of the noise, in this case generated as a random superposition of 200 pulses. The 
signal, after filtering by oblique projections along W"*", is shown in the right graph (it also 
coincides with the theoretical one). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




Figure 2: Top left graph: motion of the harmonic oscillators system as a function of time (in 
seconds) corrupted by 50 random pulses. The graph on the right depicts the signal after filtering 
the noise by sequential oblique projection. The bottom figures have the same description but 
the noise corresponds to 200 random pulses. 

Let us point out that an alternative way of splitting the signal would entail fitting both 
the signal and noise models, with the consequent increment in the dimension of the problem 
of determining the corresponding unknown parameters. In the example of this section, for 
instance, 400 more parameters (coefficients of the noise model) would be involved. 
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5 Conclusions 



Recursive equations for updating/downdating oblique projectors have been proposed. The 
updating strategy can be regarded as a generahzed Gram-Schmidt hke procedure for generating 
a sequence giving rise to oblique projectors along a fixed subspace. The downdating strategy 
modifies such sequence to account for the removal of some elements. The equations are of 
the same nature as those for producing orthogonal projectors, but involve different vectors. 
Orthogonal projectors arise within this framework as a particular case. 

The proposed technique has been applied to the problem of discriminating signals produced 
by different phenomena. The applications that have been considered assume that the signal 
model is determined by physical considerations and only the order of the model is to be adjusted. 
The task of setting the right order model is facihtated by the recursive nature of the proposed 
equations. In the two examples considered here the signal splitting is not very sensitive to the 
order of the signal model. However, an important difference between the two examples is the 
following: while in the first example an excessive overestimation of order model (maximum 
possible number of layers present in a crystal) does not prevent the extraction of the diffraction 
peak from the background, an excessive overestimation of the order of the model in the second 
example (maximum possible frequency of an oscillator) may produce the failure to separate 
the signal from the impulsive noise. The reason being that for very high frequencies the angle 
between the signal subspace and the noise subspace becomes very small, which generates an ill 
posed problem. 

The recursive feature of the proposed equations turns out to be even more important in those 
situations in which the signals splitting is achieved by stepwise selection of each component of 
the signal model. This is the subject of a recent work [33], where the present approach is shown 
to be of significant assistance. 
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Appendix A: Proof of Property [T] 

Let us recall that if vectors {vi}^^^ are linearly independent, all the duals {■u*^"'"^}^^/ are gen- 
erated by the recursive equation ( fT3l) . We need to show that such vectors satisfy: 

{vm,u'l^^) = 6m,i, m,i = 1, . . . ,k + 1. 

Proof. For A; = the relation holds because u\ = p^fp W'^iW'^ — {'^i^'^i) ~ {'Vi, Pw^-Vi). 
Therefore {v,,u\) = ;-^--^;-;-^-^wx^i) = i. 
Assuming that for /c + 1 = / it is true that 

{vm,u^i) = 5m,i, m,i = l,...,l 

we will prove that 

{vm,u^j^'^) = m,i = 1, ...,/ + 1. 

For this we need to consider four different situations with regard to the indices. 

I) m = 1, . . . , / and i = 1, . . . ,1. 

In this case {vm,(li+i) = (C.f. f|T4|) ). Hence, from the recursive equation f|T3|) . we have 

{Vm, u[^^) = {Vm, m') + = Sm,i- 

II) m = I + 1 and i = 1, . . . ,1 

Now 



so that, since = ||gfc+i|P (C.f. (lISl) 



III) m = / + 1 and i = l + l. 
This implies 



IV) m = 1, . . . , / and i = l + l. 
In this case 





2 - {vi+i,qi+i) 




llft+il 


2 



(^^m, Vi) = — ^ = 1. 



From I) II) III) and VI) we conclude that 

{Vm,u[^^) = 5m,h m,i = 1, ...,/ + 1, 



which proves that the vectors generated through (11311 are biorthogonal to vectors Vm^'m = 
l,...,fc + l. □ 
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